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R.  P.  Wei  and  G.  W.  Simona 
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INTRODUCTION 

Recent  fracture  mechanics  and  surface  chemistry  studies  of  en¬ 
vironment  assisted  crack  growth  In  gaseous  environments  have 
shown  that  crack  growth  may  be  controlled  In  eons  systems  by  the 
rate  of  surface  reactions  and  In  others  by  the  rate  of  transport  of 
the  aggressive  environment  to  the  crack  tip.  Based  on  considerations 
of  surface  reactions  and  gas  transport,  a  model  for  surface  reaction 
and  transport  controlled  fatigue  crack  growth  In  single  component 
gaseous  environments  was  proposed  and  experimentally  verified. *”* 

This  model  Is  able  to  account  for  the  Influences  of  gas  pressure  and 
cyclic  load  frequency  on  the  rate  of  fatigue  crack  growth.  In  prac¬ 
tice,  however,  there  is  usually  more  fh»n  one  gas  in  a  given  en¬ 
vironment,  The  various  component  gases  can  compete  for  surface  re¬ 
action  sites  and  therefore  alter  the  fatigue  crack  growth  response. 
The  effect  of  this  competition  needs  to  be  considered.  In  this  pa¬ 
per,  the  development  of  the  model  for  surface  reaction  and  transport 
controlled  fatigue  crack  growth  in  single  component  gaseous  environ¬ 
ments  is  briefly  reviewed.  Extension  of  the  same  considerations  to 
fatigue  crack  growth  in  a  binary  gas  mixture,  in  which  one  of  the 
components  acts  as  an  inhibitor,  is  described.  Quantitative  appli¬ 
cation  of  this  model  to  the  >.onslderatior.  of  oxygen  on  fatigue 
crack  growth  in  humid  air  is  discussed. 


MODELING  OF  FATIGUE  CRACK  GROWTH  IN  ONE  COMPONENT  GAS 

Modeling  of  environment  assisted  fatigue  crack  growth  in  single¬ 
component  gaseous  environments  was  based  on  the  proposition  that  the 
rate  of  crack  growth  in  an  aggressive  environment,  (da/dN)e,  is 


composed  of  the  sum  of  three  components.5*8*10 

(da/dN)e  -  (da/dN)r  +  (da/dN)cf  +  (da/dN)acc  (1) 

(da/dN)r  is  the  rate  of  fatigue  crack  growth  in  an  Inert  environment 
and,  therefore,  represents  the  contribution  of  "pure"  (mechanical) 
fatigue.  This  component  is  essentially  independent  of  frequency  at 
temperatures  where  creep  is  not  important.  (da/dN)acc  is  the  con¬ 
tribution  by  sustained-load  crack  growth  (that  is,  stress  corrosion 
cracking)  at  K  levels  above  Kiscc.  1  The  contribution  by  sustained- 
load  crack  growth,  t.e.  ,  the  (da/dN)scc  term,  has  been  examined  in 
some  detail  previously,11’12  and  appears  to  be  adequately  accounted 
for  by  the  superposition  model.11  (da/dN)cf  represents  a  cycle- 
dependent  contribution  requiring  synergistic  interaction  of  fatigue 
and  environmental  attack,  and  was  not  considered  by  Wei  and  Landes.11 
This  model,5  therefore,  was  directed  specifically  to  the  cycle- 
dependent  term,  and  was  limited  to  establishing  a  formal  framework 
for  estimating  frequency  and  pressure  dependence  of  (da/dN)cf  in 
single-component  gaseous  environments.  (Examinations  of  the  influ¬ 
ences  of  other  loading  variables  (such  as,  load  ratio)  are  in  prog¬ 
ress.) 

In  the  model,5  environmental  enhancement  of  fatigue  crack 
growth  is  assumed  to  result  from  embrittlement  by  hydrogen  that  is 
produced  by  the  reaction  of  hydrogenous  gases  ( e.g . ,  water  vapor) 
with  the  freshly  produced  fatigue  crack  surfaces.  More  specifically, 
(da/dN)cf  is  assumed  to  be  proportional  to  the  amount  of  hydrogen 
produced  by  the  surface  reaction  during  each  cycle,  which  is  pro¬ 
portional  in  turn  to  the  "effective"  crack  areai/produced  during  the 
prior  loading  cycles  and  to  the  extent  of  surface  reaction.  The 
time  available  for  reaction  is  assumed  to  be  equal  to  one-half  of 
the  fatigue  cycle.5  Transport  of  gaseous  environments  to  the- crack 
tip  is  assumed  to  be  by  Knudsen  flow13  at  low  pressures.  The  rates 
of  hydrogen  diffusion  and  embrittlement  were  assumed  to  be  much 
faster  than  those  of  gas  transport  and  surface  reaction,  and,  there¬ 
fore,  did  not  need  to  be  considered. 

The  governing  differential  equations  for  flow  and  surface  reac¬ 
tions,  and  the  relationship  for  (da/dN)c^  are  as  follows:5 


In  the  original  derivation,5  it  was  assumed  that  only  the  area 
produced  by  thu  previous  cycle  of  crack  growth  contributed  to  em¬ 
brittlement.  It  is  recognized  that  a  greater  portion  of  the  crack 
surface  remained  active,  particularly  for  slow  surface  reactions 
and  at  high  frequencies,  and  can  contribute  hydrogen  to  the  em¬ 
brittlement  process.8  To  provide  consistency,  while  retaining 
the  form  of  the  original  model,  an  "effective"  crack  area  or  "ef¬ 
fective"  crack  increment  La*  is  defined  here. 
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The  terms  in  the  equations  are  as  follows:5 

F  -  8.72  x  102  B(aya/E)2  (T/M)1/2  B£  (in  m3/s)  « 

Knudsen  flow  parameter  that  depends  on  dimension  and  shape 
of  the  capillary,  molecular  weight  (H)  of  the  gas  and  tem¬ 
perature  (T) .  The  specific  form  of  this  expression  re¬ 
flects  an  attempt  to  account  for  constriction  in  ilow  by 
the  real  crack,  where  £  is  a  selected  distance  (of  the 
order  10'6  a)  from  the  crack  tip  used  in  defining  a  crack 
opening  and  8  is  an  empirical  quantity  to  be  determined 
from  the  crack  growth  data.5*” 

k  *  reaction  rate  constant;  (Pa-s)~  . 

c  -2 

Nq  *  density  of  surface  sites;  molecules  (atoms) -m  . 

p  -  pressure  of  gas  at  the  crack  tip;  Pa. 

p  «  pressure  of  gas  in  the  surrounding  environment;  *a. 

3  -23  3  -1  -1 

R  *  gas  constant  -  1.38  x  10  Pa-m  -molecules  -  K 

5  ■  area  of  "effective''  crack  surface  per  cycle  ■  a  (2BAa*), 

where  Aa*  ■  "effective"  increment  per  cycle,  B  ■  specimen 
thickness,  and  a  *  empirical  constant  for  surface  rough¬ 
ness  and  crack  geometry.5 

T  ■  absolute  temperature. 

V  -  control  volume  at  the  crack  tip,  i.e. ,  volume  associated 
with  the  distance  &. 

6  *  fractional  surface  coverage  or  extent  of  reaction  of 

surface  per  unit  area. 

From  eqn.  (2),  it  can  be  seen  that  the  rate  of  change  of  pressure  at 
the  crack  tip  depends  on  the  decrease  in  pressure  produced  by  reac¬ 
tion  of  the  environment  with  the  active  ("effective")  crack  surface 
and  on  the  increase  in  pressure  from  the  influx  of  gas  from  the  ex¬ 
ternal  environment.  The  form  f(9)  ■  1  -  9  in  eqn.  (3)  incorporates 
the  assumption  that  the  surface  reaction  is  first-order  in  relation 
to  available  surface  sites.4-7 

Solutions  for  eqns.  (2)  and  (3),  with  f(9)  *1-9,  were  ob¬ 
tained  for  two  limiting  cases,  for  0  <  9  <  1,  and  were  used  in 


conjunction  with  eqn.  (4)  to  estimate  the  influences  -of  key  environ¬ 
mental  and  loading  variables  on  the  cycle-dependent  component  of  fa¬ 
tigue  crack  growth  rate,  t.e, ,  on  (da/dN)c^.5 

Case  I:  Transport  controlled. 
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Case  11:  Surface  reaction  controlled. 
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In  the  transport  controlled  case,  because  of  the  rapid  reactions  of 
the  environment  with  the  freshly  created  crack  surfaces  (high  kc)  and 
the  limited  rate  of  supply  of  the  environment -to  the  crack  tip,  sig¬ 
nificant  attenuation  of  gas  pressure  takes  place  at  the  crack  tip. 

The  extent  of  surface  reaction  (8)  during  one  cycle  is  controlled  by 
the  rate  of  transport  of  the  aggressive  environment  to  the  crack  tip, 
and  thus  varies  linearly  with  time  (see  eqn.  (5)).  i**or  the  surface 
reaction  controlled  case,  the  reaction  rates  are  sufficiently  slow 
so  that  the  gas  pressure  at  the  crack  tip  is  essentially  equal  to  the 
external  pressure.  The  extent  of  reaction,  for  f(0)  -  1  -  9,  becomes 
an  exponential  function  of  time,  eqn.  (6). 

As  a  modification  co  the  original  model,5*8  the  "effective" 
crack  increment  (Aa*)  is  now  chosen,  alc«g  with  Z  in  the  expressf-m 
of  F, to  be  equal  to  the  growth  Increment  per  cycle  at  saturation 
(•£.«.,  for  9  =  1,  Aa*  -  l  ■  (da/dN)e  8»1).  By  taking  t  equal  to  t/2 
or  l/2f,5*e  the  following  expressions  can  be  obtained  from  eqns.  (4), 
(5)  and  (6): 


Case  1:  Transport  controlled. 
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Case  II;  Surface  reaction  controlled. 


(d./dN),, 

<d.7di»  -  1  -  “P(-Vo/2f> 

*1* 


(d*/dS)  . 

TdiTSfj  —  ■  1  •'  H'<-Vo,zf) 
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(11) 


Subscript  s  is  used  to  denote  the  corresponding  values  at  saturation. 
The  modification  affects  only  the  form  of  the  expressions  for  the 
surface  reaction  controlled  case,  and  provides  a  physically  more  ac¬ 
ceptable  model.  In  their  present  form,  the  model  appears  to  be  In 
excellent  agreement  with  the  experimental  observations^/ . s_9  ’ 1  k  An 
example  of  this  correlation  is  illustrated  in  Fig.  1  for  transport 
controlled  crack  growth  in  a  2219-T851  aluminum  alloy.6 


mUSUK/UFMOUINCt  <*-,) 


Fig.  1.  Influence  of  water  vapor  pressure  (or  pressure/ 2xfrequency) 
on  fatigue  crack  growth  rates  in  2219-T851  aluminum  alloy 
at  room  temperature.  Solid  lines  represent  model  predic¬ 
tions  . 5  * 6 


27  2 

-'There  is  a  discrepancy  of  10  in  exposure  (p0/2f)  between  the 
model  predictions  and  the  experimental  data  on  AISI  4340  steel,9 
suggesting  that  capillary  condensation  of  water  vapor  may  have 
taken  place.  With  capillary  condensation,  the  effective  exposure 
can  be  much  higher  than  that  given  by  p0/2f. 


Fig.  2.  Schematic  illustration  and  comparison  of  gas  transport  and 
surface  reaction  controlled  fatigue  crack  growth. 


The  engineering  significance  and  implications  of  the  model  may 
be  discussed  in  relation  to  eqns.  (7)  and  (11),  and  considered 
through  the  schematic  illustration  shown  in  Fig.  2.  The  model  pro¬ 
vides  a  formalism  for  "predicting”  environmentally  assisted  fatigue 
crack  growth  response  in  relation  to  gas  pressure  and  cyclic  load 
frequency.  Specifically,  eqn.  (7)  gives,  for  the  transport  con¬ 
trolled  case,  the  value  of  p0/2f  at  which  "saturation"  in  environ¬ 
mental  effect  can  be  expected  in  terms  of  the  properties  of  the  al¬ 
loy,  molecular  weight  of  the  aggressive  gas,  and  temperature.  Fa¬ 
tigue  crack  growth  response  is  given,  for  this  case,  as  an  explicit 
function  of  p0/2f  below  saturation  by  eqns.  (8)  and  (9).  For  the 
surface  reaction  controlled  case,  the  dependence  of  fatigue  crack 
growth  response  on  p0/2f  is  given  by  eqns.  (10)  and  (11).  It  '*nould 
be  noted  that  eqns.  (8)  to  (11)  deal  only  with  relative  rates.  The 
actual  growth  rates  depend  on  the  Interactions  of  the  embrittling 
specie  (hydrogen)  with  the  alloy,  which  are  not  adequately  under¬ 
stood  at  this  time. 

As  further  illustration,  two  hypothetical  cases  —  one  repre¬ 
senting  gas  transport  and  the  other  surface  reaction  control  —  are 
shown  in  Fig.  2.  Crack  growth  response  curves  are  given  in  terms  of 
the  ratio  (da/dN)e/(da/dN)r ,  as  functions  of  p0/2f,  which  can  be 
readily  obtained  from  eqns.  (1),  (8)  and  (11).  In  addition,  fre¬ 
quency  scales  are  shown  for  an  external  pressure  (pQ)  of  1  kPa.  For 
this  illustration,  (da/dN) e>s/(da/dN)r  ■  50,  and  kc  »  2  x  10^  and 
2  x  10"®  (Pa-s)"*  are  used.  Although  the  curves  differ  somewhat  in 
detail  from  experimental  data,  they  do  serve  to  illurtrate  fatigue 
crack  growth  response  for  a  high  strength  steel  exposed  to  hydrogen 


sulfide  (transport  controlled)  and  water  vapor  (reaction  controlled). 
It  can  be  seen  that  the  range  of  exposures  (p0/2f)  or  frequencies 
over  which  apparent  fraquency  dependence  may  be  observed  differ  by 
about  6  orders  of  magnitude  for  the  two  cases.  The  apparent  inde¬ 
pendence  of  frequency  or  pressure,  or  p0/2f,  does  not  imply  absence 
of  environmental  effect.  Neither  does  testing  at  high  frequencies  in 
itself  ensure  absence  of  environmental  effects.  Most  importantly, 
each  material-environment  system  must  be  considered  individually. 
Generalizations  and  extrapolations  on  the  basis  of  limited  data  can 
be  very  misleading  and  should  be  avoided. 


MODELING  OF  FATIGUE  CRACK  GROWTH  IN  BINARY  GAS  MIXTURE 


Experimental  evidence  obtained  to  date1-4' 6-9  clearly  indicates 
that  the  rate  of  environment  assisted  crack  growth  is  dependent  on 
the  rate  of  hydrogen  production  at  the  crack  surfaces.  The  model  for 
fatigue  crack  growth  developed  earlier5*  -  may  be  readily  extended-  to- 
the  consideration  of  crack  growth  in  gas  mixtures.  For  simplicity, 
the  case  of  a  binary  mixture  is  considered.  One  of  the  component 
gases  is  taken  to  be  an  inhibitor  (that  is,  a  gas  which  will  react 
with  the  clean  metal  surface  but  will  not  enhance  crack  growth) .  As 
such,  the  results  may  be  used,  for  example,  to  examine  the  influence 
of  oxygen  (an  inhibitor)  on  fatigue  crack  growth  in  humid  air,  where 
water  vapor  is  the  aggressive  component.14-16 


It  is  assumed  that  (i)  both  gases  are  strongly  adsorbed  on  the 
clean  metal  surface,  (ii)  chemical  adsorption  of  either  gas  ac  a 
given  surface  site  precludes  further  adsorption  at  that  site,  (ill) 
the  ratio  of  partial  pressures  of  the  gases  at  the  crack  tip  is  es¬ 
sentially  the  same  as  that  of  the  surrounding  environment,  and  (iv) 
no  capillary  condensation  of  either  gas  occurs  at  the  crack  tip.  In 
line  with  the  model  of  Weir  et  al.,5  the  cycle-dependent  component 
of  crack  growth  rate  in  the  gas  mixture,  (da/dN)cf>m,  is  assumed  to 
be  proportional  to  the  extent  of  surface  reaction  with  the  aggres¬ 
sive  gas  during  one  loading  cycle,  or  to  6a.  If  one  now  assumes, 
for  simplicity,  that  the  kinetics  of  surface  reaction  for  both  gases 
are  first  order  with  respect  to  pressure  and  to  available  surface 
sites,  then  the  rate  equations  for  the  surface  reactions  may  be 
written  as  follows: 


d6 

-d T  "  Va(1  -  0) 

d6i 

IT  "  Vi(1  -  e) 


(12) 

(13) 


where  the  subscripts  a  and  i  denote  the  aggressive  and  inhibitor 
gases  respectively.  The  quantities  ka,  pa,  and  pf  are  the 


reaction  rata  conatants  and  partial  pressures  of  the  ‘gases  at  the 
crack  tip,  respectively.  The  coverages  8*  and  8i  denote  the  fraction 
of  surface  that  has  reacted  with  the  aggressive  and  inhibitor  gases, 
respectively.  Coverage  for  both  gases  is  denoted  by  6,  where  6  ■ 

8a  +  8i  and  0  <  8  <  1. 


Eqns.  (12)  and  (13)  may  be  solved  straightforwerdly  to  obtain 
the  extent  of  reaction  of  a  fresh  surface  with  each  gas  after  being 
exposed  to  the  gas  mixture  for  a  time  t. 


8 

a 


k  P 


V/’Vi'1 '  “p!'(k*p.+  kipi>t!> 


(14) 


8i- 


klpi 


V,+  kiPi 


-{1  -  exp[-(kgpg+  U±p±)  tl  > 


(15) 


It  follows  that  the  fraction  of  surface  that  has  reacted  with  the 
aggressive  gas  (8a)  in  relation  to  the  total  surface  coverage  (6) 
may  be  given  by  eqn.  (16). 
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If  the  total  pressure  of  the  gas  mixture  at  the  crack  tip  (pm  - 
Pa  +  Pi)  and  cyclic-load  frequency,  or  Pm/2f,  are  such  to  produce 
"complete"  surface  reaction  or  "saturation"  coverage  (8  »  1)  during 
one  cycle,  then  the  fraction  of  the  crack  surface  that  has  reacted 
with  the  aggressive  gas  is  given  by  eqn.  (17) . 
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Since  (da/dN) *  8*Aa*  (see  eqn.  (4)),  the  cycle-dependent 
component  of  fatigue  crack  growth  rate  in  the  mixed  gases, 

(da/dN) cf  m,  can  be  obtained  from  the  cycle-dependent  rate  in  a  pure 
gas  at  saturation,  (da/dN)  _  ,  and  eqn.  (17) . 
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1) 


The  subscript  o  Is  used  to  denote  pressures  in  the  external  environ¬ 
ment,  and  r***lects  the  assumption  that  the  partial  pressure  ratio  at 
the  crack  tip  is  the  same  as  that  in  the  external  environment.  It 
should  be  noted  that  only  the  competition  between  two  gases  has  been 
modeled  here,  and  eqn.  (18)  is  expected  to  apply  equally  well  to 
both  transport  and  surface  reaction  controlled  crack  growth.  Sche¬ 
matic  illustration  of  eqn.  (18)  is  shown  in  Fig.  3  as  (da/dN)e  m/ 
(da/dN)r  for  different  ratios  of  (da/dN)ejS/(da/dN)r;  the  ratios 
reflect  different  degrees  of  severity  of  environmental  effects  in 
the  pure  gas.  A  ratio  of  ki/ka  of  0.14,  which  approximates  the  ratio 
between  the  reaction  rate  constants  for  oxygen  and  water  vapor  over 
a  2219-T851  aluminum  alloy,6'17  is  used.  The  upper  scale  for  pa  in 
Fig.  3  is  computed  for  pi  *  152  torr  (or  20.2  kPa) ,  and  therefore 
corresponds  to  the  partial  pressure  of  water  vapor  in  humid  (moist) 
air. 


Assuming  that  the  surface  reaction  data  on  2219-T851  aluminum 
alloys6*17  may  be  used  for  other  high-strength  aluminum  alloys,  com¬ 
parisons  of  model  predictions  with  data  reported  by  Hartman  et  al.15 
and  Feeney,  McMillan  and  Wei16  on  7075-T6  aluminum  alloy,  and  Brad¬ 
shaw  and  Wheelei  1 4  on  DTD  5070A  aluminum  alloy  are  made  and  are 
shown  in  Figs.  4  and  5,  respectively.  For  the  7075-T6  aluminum  al¬ 
loy,  data  in  water  and  salt  water16  are  assumed  to  represent  the 
"saturation"  level  m  pure  water  vapor  (8  -  l).6  Excellent  agree¬ 
ment  is  observed  at  the  higher  K  levels  and  higher  crack  growth 
rates.  At  the  lower  K  levels  and  lower  crack  growth  rates,  the  data 
tend  to  follow  those  of  pure  water  vapor,  and  may  reflect  partial 
capillary  condensation  during  a  portion  of  each  loading-unloading 


Fig.  3.  Schematic  illustration  of  the  influence  of  inhibitor  gas 
on  fatigue  crack  growth  in  a  binary  gas  mixture.  (In  the 
upper  scale,  pa  corresponds  to  pi  ■  152  torr  or  20.2  kPa.) 


WATER  VAPOR  PRESSURE  (torr) 


Fig.  4.  Influence  of  water  vapor  pressure  on  fatigue  crack  growth 
in  7075-T6  aluminum  alloy  in  humid  air.15’16  Dashed  and 
solid  lines  represent  model  predictions  from  eqns.  (7), 
(8)  and  (18). 


PRESSURE/ 2  *  FREQUENCY  (Po-*) 


Fig.  5.  Influence  of  water  vapor  pressure  and  frequency  on  fatigue 
crack  growth  in  DTD  5070A  aluminum  alloy.14  Solid  lines 
represent  model  predictions  for  water  vapor,  and  dashed 
lines  for  humid  air. 


cycle  or  some  other  factor.  Further  studies  are  needed  to  clarify 
this  discrepancy. 


DISCUSSION 


On  the  basis  of  available  experimental  results ,  it  is  seen  that 
the  models  reasonably  account  for  environmentally  assisted  fatigue 
crack  growth  response  in  single-component  gaseous  environments  and 
in  binary  gas  mixtures,  in  which  one  of  the  components  acts  as  an 
inhibitor.  A  similar  consideration  of  environmentally  assisted 
crack  growth  in  binary  gas  mixtures,  under  sustained  loading,  has 
also  been  made  and  provided  good  agreement  with  experimental  obser¬ 
vations.18  Thus,  the  models  provide  a  formalism  for  the  planning  of 
experiments,  and  for  the  interpretation  and  extrapolation  of  data 
for  engineering  applications. 

In  deriving  these  models  and  in  their  application,  only  simple 
surface  reactions  have  been  considered.  Many  reactions  that  are  of 
Interest  (for  example,  reactions  of  steels  with  water  vapor  or  hy¬ 
drogen  sulfide),  however,  are  more  complex.  Reactions  of  different 
components  to  form  new  complexes  are  also  possible.  Incorporation 
of  these  complexities  into  future  modeling  efforts  must  be  consider¬ 
ed.  Extension  of  these  ideas  to  the  consideration  of  fatigue  crack 
growth  in  mixtures  of  H2S  and  CO  has  been  made.19'20  Preliminary 
results  Indicate  that  modeling  of  the  more  complex  environments 
would  be  feasible. 


SUMMARY 

Based  on  considerations  of  surface  reactions  and  gas  transport, 
models  of  environmentally  assisted  fatigue  crack  growth  in  single¬ 
component  gaseous  environments  and  in  binary  gas  mixtures  (in  which 
one  of  the  components  acts  as  an  inhibitor)  have  been  developed. 

Good  agreement  with  experimental  data  has  been  observed.  The  models, 
therefore,  provide  a  formalized  framework  for  experimental  design, 
and  for  the  interpretation  and  extrapolation  of  data  for  engineering 
applications.  Further  research  will  be  needed  to  refine  and  extend 
this  work  to  the  consideration  of  more  complex  environments,  and  to 
incorporate  the  influences  of  metallurgical  variables. 
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